function [pi] = pi_fun(P,I,sigma,N,f)

for i=1:N
    pi(:,i) = I(1,i).*P.^(-sigma)-f;
end

end
